A shrinking Compact Symmetric Object: J11584+2450? 
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ABSTRACT 

We present multi-frequency multi-epoch Very Long Baseline Array (VLBA) 
observations of J11584+2450. These observations clearly show this source, previ- 
ously classified as a core-jet, to be a compact symmetric object (CSO). Compar- 
isons between these new data and data taken over the last 9 years shows the edge 
brightened hot spots retreating towards the core (and slightly to the west) at ap- 
proximately 0.3c. Whether this motion is strictly apparent or actually physical 
in nature is discussed, as well as possible explanations, and what implications a 
physical contraction of J11584+2450 would have for current CSO models. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: individual 
(J11584+2450) — galaxies: nuclei — galaxies: jets — radio continuum: galaxies 



1. Introduction 

Compact symmetric objects (CSOs) are now a well established class of radio sources 
loosely defined as sources with emission on both sides of the core (which itself is not always 
detected) on a size scale of 1 kpc or less (j Wilkinson et al.lll994j ). The generally accepted 
explanation for the small size of the se objects is that they are young radio sources which 
could grow into larger FR II objects (IReadhead et al.lll996l ; |OJDeaJll9j)8|). Alternately, it has 
been proposed that the s mall size of these structures is due to their growth being frustrated 
by a dense environment (Ivan Breugel et al.l 1 19841 : lO'Deal Il998l ) . 
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Due to their rapid growth, age estimates for the emission from these objects can be 
obtained kinematically (lOwsianik fc Conwaylll998l ; iTaylor et al.ll2000f ) yielding ages ranging 



from tens to thousands of years. Less acc urate spectroscopic models ( IReadhead et al.lll996 



Murgia et al.lll999l ; iGugliucci et al.ll2005l ) place CSOs at a few thousands of years old. This 
generally supports the theory that these are young AGN in the early stages of evolution. 
However, sinc e the age distribution of currently known CSOs is heavily weighted on the 
younger side (IGugliucci et al.l 120051 ) . this indicates the evolution might not be straightfor- 
ward. This distribution should not be taken as definitive though, since Gugliucci et al. 
concede it may be influenced by selection effects. One theory suggests CSOs could be gener- 
ally short-lived objects with only a small fraction of them surviving t o become larger-scale 
objec ts while the remaining galaxies become permanently radio quiet (lOwsianik &: Conway 
19981 ). In a competing theory there exists a cyclic process where unsucce ssful CSOs have 
multiple opportunities to grow into larger objects (jO'Dea fc Bauml 119971 ). Alternatively, 
the current distribution could be an artifact of the small statistical sample from which it is 
derived combined with selection effects. 

The picture presented in the above models might be overly simplistic. For example, 
all of thes e models predict c ontinuous radial expansion of the lobes, but CSOs such as 
1031+567 (ITaylor et al.ll2000l ) have been observed with non-radial motion. Here we present 
observations of J11584+2450 (P KS 1155+251 SP SS J115825.79+245018.0), a galaxy with 
a redshift of 0.20160 ± 0.00040 (IZensus et al.ll2002l ) that appears to be contracting towards 
its core on both sides. 



Throughout this discussion, we assume H =73 km s 1 Mpc 
so 1 mas = 3.213 pc. 
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2. Observations and Data Reduction 

Multi-frequency observations of J11584+2450 were performed on September 19, 2006 
with the VLBA. A summary of these and other observations referred to in this paper is 
presented in Table [TJ These observations consisted of four 8 MHz wide IFs in the C, X 
and U bands with full polarization centered at: 4605.5 , 4675.5, 4990.5, 5091.5, 8106.0, 
8176.0, 8491.0, 8590.0, 14902.5, 14910.5, 15356.5 and 15364.5 MHz at an aggregate bit 
rate of 256 Mbps to maximize (u,v) coverage and sensitivity. When the data in each band 
were combined, the three central frequencies were: 4844.7, 8344.7, and 15137.5 MHz. The 
integrations were performed in blocks (~2 minutes for 5 and 8 GHz, ~7.5 minutes for 15 
GHz) and these blocks were spread out over a 9.5 hour period to maximize (u,v) coverage 
of the source. 
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Most of th e calibration and initial imagi ng; of the new data were carried out by au- 



tomated AIPS (jGreisenl 120031 ) and DIFMAP (Shepherd 1997 ) scripts similar to those used 



in reducing the VIPS 5 GHz survey data (IHelmboldt et al.l 



20071 : iTavlor et al.l 120050 . To 



summarize, flaggi ng of bad data and calib ration were performed using the VLBA data cal- 
ibration pipeline ( Siouwerman et all 2005 ). while imaging was performed using DIFMAP 
scripts described in ITavlor et al.l ( 2005 ). Final imaging was performed manually using the 
DIFMAP program, with beam sizes of 1.906 x 3.16 in position angle -6.89°, 1.195 x 1.788 
in position angle -3.58° and 0.6876 x 0.9794 in position angle -2.717° for 5, 8 and 15 GHz 
respectively. 



3. Results 
3.1. Images 



Kellermann et al.l (120041 ) observed J11584+2450 (B1155+251) as part of the VLBA 2cm 
Survey. Since these observations were only at 15 GHz, they typically identified the brightest 
component in an image to be the core, and consequently classified this source as a core-jet 
with the core being the southern, bright component. 

Figure [1] shows the 5, 8, and 15 GHz VLBA images made from the September 2006 
observations of J11584+2450. The 15 GHz map shows the clearest structure so is used to 
label components of the source. The 15 GHz image shows a compact unresolved component 
(C) with resolved emission both to the north (Nl) and the south (SI). The southern emission 
then seems to have another component that expands out towards the west (W2). There also 
exists some weak emission on the western edge of the image (Wl; with a peak flux density 
of 0.55 mJy/beam). SI is the brightest component in the image (52.3 mJy/beam peak), 
and is what was previously identified as the core. In the 8 GHz image the edges of C, Nl, 
SI, and W2 become indistinguishable, but these components can still be identified by local 
peaks within the image. Interestingly, an eastern spur develops from the southern edge of 
Nl, extending in the opposite direction from the majority of the diffuse emission. Overall 
the emission appears more extended, and Wl has a more significant detection. The 5 GHz 
map further smears the interior components together until only Nl and SI are clearly visible 
as local maxima of the map. The spur mentioned above becomes brighter, and the western 
emission stretches out farther towards a very well detected Wl (3.83 mJy/beam peak flux 
density) . There is emission from the eastern spur towards the south in this image (hereafter 
referred to as the 5 GHz southeastern clump), which has an integrated flux density of 3.78 
mJy 
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The geometry between Nl, C, and SI was measured using the 2006 15 GHz data, since 
those components are most distinguishable. The axial ratio Nl/Sl = 1.59, and the angle 
subtended between the arms, Nl-C-Sl, is 166.9°. 

We used the VLA in the D- Configuration to investigate what appeared to be an exten- 
sion of the western emission to kilo-parsec scale from the NRAO VLA Sky Survey (NVSS; 
Condon et al. 1998), but found no indication of any western emission from Jl 1584+2450 
(Fig. [2]) and the previous extension to be a result of the higher RMS (0.45 mJy/beam) of 
the NVSS compared to our image with RMS = 0.085 mJy/beam. 

Additionally, we acquired a visual image (Fig. [3]) from the Sloan Digital Sky Survey 
Data Release 5 (SDSS DR5; Adelman-McCarthy et al. 2007 ). This shows the source to have 
a galaxy 4 arcseconds to the south east, which is uncatalogued outside of the SDSS. Two 
different algorithms have been used to determine a mulit-color photometric redshift for this 
object from the SDSS imag e. The first algorit hm utilized the template fitting method and 
yields z = 0.0008 ± 0.0226 (ICsabai et al.l 120031 ) , while the second algorithm used a Neural 
Network method and yields z = 0.065 ± 0.123 (jOyaizu et al.l 120071 ). placing an upper limit 
of z 0.19 which is comparable to J11584+2450's redshift. 



3.2. Spectral Index Distribution 

The 2006 VLBA images at 8 and 15 GHz were matched in resolution in order to obtain 
a spectral index distribution across the source that was overlaid onto a 5 GHz image to show 
overall source structure (Figure H]). This distribution clearly shows a compact flat-spectrum 
component (a ~ —0.276, where F u oc u a ) situated between two steeper spectrum lobes 
(a ~ —1.03). More steep spectrum emission is found to the west of these lobes (a ~ —1.36 
to —1.54), where it then fades below the detection threshold at high frequencies. Looking 
at the 5 GHz southeastern clump and using its peak, the spectral index would have to be 
steeper than —4.35 for the emission to fall below the RMS of the 8 GHz image. Alternatively, 
the absence of this feature could be due to having fewer short spacings at 8 GHz, or it could 
merely be an imaging artifact at 5 GHz. 
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4. Discussion 



4.1. Reclassification of J11584+2450 



The compact flat spectrum c omponent seen in Figu re S] is compatible with emission from 
the nucleus, or core, of a galaxy (IBegelman et al.lll984l ). The steeper spectrum components 
extending north and south from the core are accordant with the spectral signature of jets or 
hot spots. This overall structure is clearly consistent with J11584+2450 being a CSO. 



4.2. Component Motions 

To characterize intrinsic motions within the source, a multi-component elliptical gaus- 
sian model was made of the 1999 data and then applied to the 1995, 2001 and 2006 15 GHz 
data varying the flux and position parameters of each gaussian component (Table [2]). The 
extended dual-lobed structure in the 1995 data (but notably missing from the 1999, 2001 
and 2006 images) of what is now considered the core was modeled using a single component, 
since it is likely this extension is only an artifact. 

Each of the four IFs of the 2006 data were then individually modelled and and com- 
pared to each other to determine the systematic error (<r sys ) associated with modelling each 
component. The total position error (a tot) wa s then calculated for individual components 
using of ot = cr^ tat + cr^ s , where a sta t is the expected statistical error associated with mod- 
eling gaussian components in two dimensional polar coordinates (Table [3]) (adapted from 
1-dimension case per Fomalont 1999). Since all model positions are referenced to C, a sys for 
the core is accounted for by the other component uncertainties. Since the calibrated data 
sets from the VLBA 2cm Survery are each averaged to one frequency, the a sys attained from 
the 2006 data was applied to them as well. Similarly, a multi-component gaussian model 
was made of the 2000 data and then applied to the 2006 8 GHz data varying the flux and 
position parameters of each gaussian component (Table 0J. a sys values were obtained for 
each component using the four IFs as above in the 2006 data, and using the four IFs closest 
to those same frequencies in the 2000 data (Table [5]). 

These models were used to calculate component velocities relative to the core for each 
band, which are plotted in Fig. [5] with the tail of each vector located at the earliest data 
position in the band. The components representing the hot spots at the working surface of 
the jets which were modeled (Nl and SI) appear to have contracted towards the core as well 
as traveled westward, and this motion is consistent between the two bands. Additionally, the 
W2 component moves towards the northwest in both bands. Performing a least-squares fit 
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to solve for the velocities in each band separately, and then using these independent values 
to reduce the error yields a radial contraction velocity (normalized to the speed of light) of 
0.42 ±0.03 c for SI and 0.20 ±0.07 c for Nl. Overlaying the contour maps of different epochs 
(see Fig. [6] for one example) is also supportive of contraction, since the 1995 contours are 
interior to the 2006 contours suggesting that the component motion is not an artifact of the 
modeling. 

This motion was not detected by the VLBA 2 cm Survey over the six year interval 
between 1995 and 2001 since the velocities i nvolved (~ 0.03 mas /yea r) are well within most 



of their stated velocity errors for this source (IKellermann et al.ll2004l ). The total flux density 



of the source has been decreasing steadily since at least 1995, the models show that this 
drop can be attributed to SI decreasing steadily (40% decrease at 15 GHz over this 11 year 
period), while the other components exhibit small fluctuations. 



4.3. Apparent Motion Interpretation 

Since actual contraction of the source towards the core is something that has not been 
previously observed, we first examine reasons behind an effect that would merely cause 
apparent motion in the system. One possible explanation for seeing the contraction of 
this source is that if hotspots are advancing out away from the core and expanding and 
younger hotspots are brightening due to interactions at the end of the jet, then the models 
might not be fit to the same components. The largest problems with this hypothesis are 
its lack of explanation of both the western emission and the western component to the 
hotspot velocities, which means these properties require a separate unrelated explanation if 
the contraction is to be explained by hot spot dimming and advance. 



4.4. Physical Motion Interpretations 



Leaving open the possibility that the data represent physical motions in the system, we 
include discussion along those lines. One interpretation is that we are viewing a projection 
effect caused by rotation of the source. While solid body rotation is an unphysical scenario, 
it gives us an idea about what fluid rotation would look like for this system so we consider 
it as a first approximation of rotational motion. The angular velocity of the rotation is 
dependent upon the initial orientation. Assuming the axis of rotation lies in the plane of 
the sky, and the inclination angle is less than 45° since larger values would yield Doppler 
boosting, the jets would have a rotational period between just 260 and 1880 years. 
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Since the system is a fluid and not a rigid body it would actually have differential rota- 
tion. The core would therefore be spinning even faster than the observed jet components and 
we would expect the core to appear highly variable since the base of the jet would frequently 
be pointed towards us. A second physical interpretation of the motion is precession. The 
westward component to the velocities of both jets strongly argues against precession. If the 
jets were precessing and thus appearing shorter, they would move in opposite directions (i.e. 
one moves east while the other moves west) as well as inward. Additionally, these hypotheses 
fail to explain the older western emission. 

Another physical elucidation of the motion is that there exists some reason for the 
pressure of the environment to increase, then this could leave the jets under-pressured and 
lead to contraction of J11584+2450. Such a pressure change could result from a relative 
motion between a clumpy environment and the CSO. The 1.59 Nl/Sl axial ratio is also 
indicative of the jets encountering a dense environment. If the departure of the ratio from 1 
was due to Doppler boosting then the brighter hotspot would be further from the core (SI 
in the case of J11584+2450). This would yield a ratio smaller than 1, therefore the ratio is 
likely due to the jet running into difficulty as it tunnels through the environment causing it 
to be both shorter and brighter. Both the angle Nl-C-Sl and the western deviation of the 
emission is consistent with the source moving eastward and being influenced by ram pressure 
similar to what has been observed in wide-angle tailed radio sources (e.g. 3C465, Hardcastle 
et al. 2005; Sakelliou & Merrifield 2000), but on a smaller spacial scale. Relative motion 
could exist between J11584+2450 and its host galaxy causing the interstellar medium to 
produce ram pressure against the radio jet. If the companion galaxy to the south east (Fig. 
[3]) is at a similar redshift and these two galaxies are members of a cluster then Jl 1584+2450 
might be moving towards the center of the gravitational potential. However, this is highly 
speculative and more observations are needed to determine the cluster environment. 



5. Conclusions 

After analyzing multifrequency (5, 8, and 15 GHz) VLBA data from radio source 
J11584+2450 we reclassify it as a CSO. Fitting the data with multi-component gaussian 
models and overlaying images of different epochs on each other not only show that this 
source is not growing at the usual rate of ~ 0.1 to 0.3 c, but each jet is apparently shrinking 
in size at ~ 0.3 c and each is additionally moving westward at ~ 0.2 c. Co nfirmation of 



other CSOs having either recessive behavior or non-radial motion like 1031+567 (ITaylor et al. 



20001 ) would mean the current models need to be modified to allow for possible non-linear 



growth periods during the evolution of AGN. The prospect of non-linear growth for CSOs 
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would further bring in to question the validity of kinematic ages. iGugliucci et al.l (120051 ) 
found 7 out of the 13 CSO they dated to be under 500 yr old and commented that the ex- 
pectation for a steady-state population of CSOs would have a uniform distribution of ages. 
While this is a small statistical sample it is also what one would expect to see if CSOs spend 
a greater fraction of time as small sources. However, in kinematic observations of ~ 10 CSOs 
we have seen contraction in just 1 source. 

Future VLBA observations of this source are planned to follow the motion of the com- 
ponents and to see whether they continue to recede towards the core, and what time-scale 
this occurs over. Lower frequency (1.4 GHz) observations should be carried out to confirm 
the existence of the 5 GHz southeastern clump. If the emission from both jets is actually 
flowing towards the west, lower frequency observations might also show Wl merging with 
the diffuse western emission, and could reveal larger scale structures. 
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Fig. 1. — VLB A observations from September 2006 of Jl 1584+2450 at frequencies of (moving 
top to bottom) 4.84, 8.34, and 15.13 GHz. Contour levels begin at 0.375 mJy/beam and 
increase by factors of 2 1 / 2 . 
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Fig. 2.— VLA observations from March 2007 centered on J11584+2450 at 1.3649 GHz in 
D configuration. The contour levels begin at 0.336 mJy/beam and increase by levels of 2. 
The source shows no significant structure at this frequency and resolution, but does show 
possible extension to the south-east. 
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Fig. 3. — This R-Band SDSS image of Jl 1584+2450 shows a second source 4 arcseconds to 
the southeast (SDSS J115826. 16+245014. 9). If these sources have the same redshift, then 
there is 13 kpc separation between them. In this image, the magnitude of Jl 1854+2450 is 
17.68 ± 0.01, while the magnitude of SDSS J115826.16+245014.9 is 21.99 ± 0.06. 
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Fig. 4. — 1. Multi- Frequency observations of a newly identified CSO (J11584+2450). 5 GHz 
contours overlaid on a 8 to 15 GHz spectral index image. Notice the flat spectrum core, as 
well as the symmetric dual-lobed structure in the source. Also, the emission abruptly bends 
to the west. This sudden path change, and the steep spectrum compact knot at the western 
edge is not clearly understood. 




Fig. 5. — Relative velocity of model components. Velocity of each gaussian model component 
is plotted (lmas = 0.2c) with the tail of each vector originating at the model components 
position at its earliest observation (1995 for 15 GHz & 2000 for 8 GHz). The model fits are 
in agreement with the contour overlay plots in showing this source to be shrinking. 
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Fig. 6. — Two epoch overlay of J11584+2450. 15 GHz data from 1995 (Red Contours) is 
plotted with the 15 GHz data from 2006 (Blue Contours). This figure shows contraction of 
the source towards the core (denoted by a cross) over time. The contour levels begin at 1.10 
mJy/beam and increase by factors of 2 1 / 2 . 



Table 1. VLA and VLBA Observations of J1158+2450. 



Freq. 


Date 


Time 


BW 


Pol. 


IFs 


Peak 


rms 


(GHz) 




(minutes) 


(MHz) 






(mjy beam ) 


(mjy beam ) 


1.3649* 


2007 Mar 08 


24.7 


100 


4 


2 


1073.7 


0.1 


4.8447+ 


2006 Sep 19 


24.9 


32 


4 


4 


193.38 


0.11 


8.3447 t 


2006 Sep 19 


26.9 


32 


4 


4 


118.38 


0.17 


8.354lt 


2000 May 06 


20.6 


64 


1 


8 


133.1 


0.2 


15.138t 


2006 Sep 19 


90.8 


32 


4 


4 


52.27 


0.11 


15.335 1 " * 


2001 Mar 04 


57.0 


56 


1 


1 


74.85 


0.31 


15.335+ * 


1999 May 21 


37.6 


56 


1 


1 


82.73 


0.33 


15.350+ * 


1995 Apr 07 


44.9 


32 


1 


1 


110.0 


0.3 



*VLA Observation 
+VLBA Observation 

■l-These data were taken as part of the VLBA 2 cm Surverv, iKellermann et all (2004) 
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Table 2. 15 GHz Gaussian Model Components* . 



Component 


Epoch 


S 


r 




6 


0-0 


o 


b/a 


4> 






(Jy) 


(mas) 


(mas) 


(°) 


(°) 


(mas) 




(°) 


C... 


1995 


0.0167 


0.000 


0.007 


0.00 


0.00 


0.283 


1.00 


-18.57 




1999 


0.0331 


0.000 


0.003 


0.00 


0.00 


0.283 


1.00 


-18.57 




2001 


0.0170 


0.000 


0.007 


0.00 


0.00 


0.283 


1.00 


-18.57 




2006 


0.0120 


0.000 


0.003 


0.00 


0.00 


0.283 


1.00 


-18.57 


SI... 


1995 


0.1267 


3.647 


0.014 


-161.28 


0.15 


0.272 


0.86 


50.59 




1999 


0.0945 


3.514 


0.014 


-159.79 


0.15 


0.272 


0.86 


50.59 




2001 


0.0798 


3.400 


0.014 


-159.03 


0.15 


0.272 


0.86 


50.59 




2006 


0.0765 


3.360 


0.014 


-160.12 


0.15 


0.272 


0.86 


50.59 


Nl... 


1995 


0.0327 


5.610 


0.052 


9.55 


1.12 


3.02 


0.345 


29.62 




1999 


0.0288 


5.801 


0.049 


8.37 


0.95 


3.02 


0.345 


29.62 




2001 


0.0257 


5.497 


0.045 


9.52 


1.78 


3.02 


0.345 


29.62 




2006 


0.0399 


5.325 


0.038 


6.82 


0.50 


3.02 


0.345 


29.62 


W2... 


1995 


0.0637 


3.226 


0.034 


-127.77 


0.46 


2.72 


0.32 


-81.22 




1999 


0.0541 


3.775 


0.033 


-115.06 


0.46 


2.72 


0.32 


-81.22 




2001 


0.0369 


3.746 


0.033 


-115.65 


2.65 


2.72 


0.32 


-81.22 




2006 


0.0630 


3.435 


0.030 


-116.83 


0.34 


2.72 


0.32 


-81.22 



*NOTE - Parameters of each Gaussian component of the model brightness distribution are 
as follows: Component, Gaussian component; Epoch, year of observation (see Table [TJ; S, flux 
density; r, oy, 8, erg, polar coordinates (and the associated errors) of the center of the compo- 
nent relative to the center of component C; a, semimajor axis; b/a, axial ratio; <3>, component 
orientation. All angles are measured from north through east. 



Table 3. 15 GHz Gaussian Model Components from Selected IFs* . 



Component 


Epoch 


ri 


T2 


f"3 


T4 


^r(sys) 


6*1 


6 2 


831 


84 


a 0(sys) 






(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(°) 


(°) 


(°) 


(°) 


(°) 


SI... 


2006 


3.368 


3.353 


3.341 


3.377 


0.014 


-159.93 


-160.10 


-160.36 


-160.13 


0.15 


Nl... 


2006 


5.276 


5.343 


5.376 


5.318 


0.037 


7.20 


6.19 


6.75 


7.19 


0.41 


W2... 


2006 


3.473 


3.430 


3.390 


3.440 


0.030 


-116.45 


-116.63 


-117.06 


-117.30 


0.34 



*NOTE - Parameters of each Gaussian component of the IF model position distribution are as follows: Component, 
Gaussian component; Epoch, year of observation (see Table [T}; ri — r^, radial positions of model components in IFs 1 
through 4 respectively; cr r ^ sys ^, the standard deviation in radial position; 81 — 84, the polar angular position of model 
components in IFs 1 through 4 respectively; tTg( sys ), the standard deviation in polar angular position. All angles are 
measured from north through east. 
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Table 4. 8 GHz Gaussian Model Components* . 



Component 


Epoch 


S 


r 




6 




a 


b/a 


* 






(Jy) 


(mas) 


(mas) 


(°) 


(°) 


(mas) 




(°) 


C... 


2000 


0.0273 


0.000 


0.001 


0.00 


0.00 


0.214 


1.00 


63.74 




2006 


0.0115 


0.000 


0.003 


0.00 


0.00 


0.214 


1.00 


63.74 


SI... 


2000 


0.1581 


3.240 


0.017 


-157.29 


0.15 


0.4875 


0.71 


-23.38 




2006 


0.1507 


3.082 


0.022 


-155.45 


0.19 


0.4875 


0.71 


-23.38 


Nl... 


2000 


0.1005 


5.597 


0.055 


6.75 


0.31 


2.195 


0.76 


-38.51 




2006 


0.1098 


5.508 


0.049 


5.27 


0.10 


2.195 


0.76 


-38.51 


W2... 


2000 


0.1505 


3.653 


0.028 


-109.43 


0.32 


1.57 


0.80 


40.08 




2006 


0.1489 


3.703 


0.011 


-105.95 


0.47 


1.57 


0.80 


40.08 



*NOTE - Parameters of each Gaussian component of the model brightness distribution are 
as follows: Component, Gaussian component; Epoch, year of observation (see Table [TJ; S, flux 
density; r, cr r , 0, cr$, polar coordinates (and the associated errors) of the center of the compo- 
nent relative to the center of component C; a, semimajor axis; b/a, axial ratio; component 
orientation. All angles are measured from north through cast. 



Table 5. 8 GHz Gaussian Model Components from Selected IFs* . 



Component 


Epoch 


n 


V2 




T4 




01 


02 


6 3 i 


04 


°e(sys) 






(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(°) 


(°) 


(°) 


(°) 


(°) 


SI... 


2000 


3.225 


3.270 


3.250 


3.237 


0.017 


-157.64 


-157.48 


-157.26 


-157.29 


0.15 




2006 


3.106 


3.079 


3.051 


3.101 


0.022 


-155.60 


-155.27 


-155.24 


-155.67 


0.19 


Nl... 


2000 


5.573 


5.675 


5.535 


5.551 


0.055 


7.31 


6.65 


7.05 


6.54 


0.31 




2006 


5.467 


5.501 


5.586 


5.467 


0.049 


5.17 


5.34 


5.19 


5.39 


0.09 


W2... 


2000 


3.607 


3.673 


3.675 


3.642 


0.028 


-109.24 


-109.98 


-110.07 


-109.45 


0.32 




2006 


3.692 


3.721 


3.696 


3.707 


0.011 


-106.46 


-105.96 


-105.23 


-106.29 


0.47 



*NOTE - Parameters of each Gaussian component of the IF model position distribution are as follows: Component, 
Gaussian component; Epoch, year of observation (see Table [TJ; r\ — t"4, radial positions of model components in IFs 1 
through 4 respectively; cr T / s y S \, the standard deviation in radial position; 0i — 04, the polar angular position of model 
components in IFs 1 through 4 respectively; og( ByB ), the standard deviation in polar angular position. All angles are 
measured from north through east. 



